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transitions caused by pulses and (ii) due to external perturbation. The scaling of errors with system 
parameters and number of qubits is explained. We use two different procedures to fight each of 
them. To suppress spurious transitions we use correcting pulses (generalized 2nk method) while to 
suppress errors due to external perturbation we use an improved QFT algorithm. As a result, the 
fidelity of quantum computation is increased by several orders of magnitude and is thus stable in a 
much wider range of physical parameters. 
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I. INTRODUCTION 

Quantum information theory pj is a rapidly evolving 
field. It uses quantum systems to process information 
and by doing so can achieve things that are not possi- 
ble with classical resources. Quantum secure commu- 
nication for instance is already commercially available. 
Quantum computation on the other hand is still far from 
being useful outside of the academic community. Two 
serious obstacles to overcome in building quantum com- 
puters are: (i) one must be able to control the evolu- 
tion in order to precisely implement quantum gates, (ii) 
one must suppress all external influences. Errors in both 
cases are caused by the perturbation of an ideal quantum 
computer, either due to the internal imperfections in the 
first case or due to coupling with the "environment" in 
the second case. In the present paper we study both kind 
of errors in QFT algorithm and try to minimize them. 

In order to be closer to the experimental situation we 
choose a concrete model of a quantum computer, namely 
the Ising quantum computer (IQC) [2|- IQC is one of 
the simplest models still having enough complexity to 
allow universal quantum computation. Quantum gates 
on this computer can be realized by the application of 
electromagnetic pulses. For the algorithm we choose to 
discuss QFT. The first reason to choose QFT is that it is 
one of the most useful quantum algorithms, giving expo- 
nential speedup over the best classical procedure known, 
and is also one of the ingredients of some other impor- 
tant algorithms, e.g. Shor's factoring algorithm 3]. The 
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second reason is that it is a complex algorithm, where 
by complex we mean it has more than 0(n) number of 
quantum gates as opposed to previously studied more 
simple algorithms where the number of gates scales only 
linearly with the size of the computer (e.g. entangle- 
ment protocol 4]). Previous study of Shor's algorithm 
m IQC || did not use recently introduced generalized 
27rfc method [|| which is the best known procedure for 
inducing transitions on IQC. It is easy to imagine that in 
most useful quantum algorithms the size of the program 
will grow faster than linearly with the number of qubits 
n and therefore it is important to see how errors accu- 
mulate in such algorithms. This importance is confirmed 
by our results showing that errors due to unwanted tran- 
sitions for QFT grow with the square of the number of 
pulses and not linearly as in algorithms with linear 0{n) 
number of gates, for a typical state. 

For QFT algorithm running on IQC we analyze er- 
rors due to spurious transitions caused by pulses, these 
we call intrinsic errors, and errors due to the coupling 
with an external "environment" , called external errors, 
modeled by a random hermitian matrix from a gaussian 
unitary ensemble (CUE) 7]. We minimize intrinsic er- 
rors by applyingsome additional pulses to correct most 
probable errors ,6] and by doing this, we are able to sup- 
press intrinsic errors by several orders of magnitude. To 
suppress external errors due to GUE perturbation we use 
previously proposed improved quantum Fourier transfor- 
mation (IQFT) Q which is more stable against GUE 
perturbations in a certain range of parameters. By using 
correlation function approach |9( we analyze in detail the 
dependence of errors on all relevant parameters and on 
the number of qubits n. By doing this we can set the lim- 
its between which parameters of IQC should lay in order 
to preserve the stability of computation. In our approach 
to decrease errors we do not use error correcting codes for 
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the following reasons: we want to remove as many errors 
as we can on the lowest possible level and second, the in- 
trinsic and external errors are not easily handled by error 
correcting codes (see Ref.[T(il] and references therein). 

The outline of the paper is as follows. In section II 
we repeat the definition of IQC and in section III we 
summarize linear response formalism which is the main 
theoretical tool for studying fidelity. In section IV we 
study intrinsic and external errors, first separately and 
then the case of both errors present at the same time. 
In the appendix we present the pulse sequences used to 
implement QFT and IQFT algorithms. 



II. ISING QUANTUM COMPUTER 

IQC consists of a 1-dimcnsional chain of n equally 
spaced identical spin 1/2 particles coupled by nearest 
neighbor Ising interaction of strength J, so that paral- 
lel spins are favored over anti-parallel ones by an energy 
difference of J (we set h = 1 throughout the paper). The 
quantum computer is operated via an external magnetic 
field having two components. The first one is a perma- 
nent magnetic field oriented in the z direction with a con- 
stant gradient which allows for the selective excitation of 
individual spins, while the second one is a sequence of 
T circular polarized fields in the x-y plane (which are 
called pulses) with different frequencies v^ m \ amplitudes 
(proportional to the Rabi frequencies r^ m '), phases <^( m ) 
and durations r'" 1 ^ for the mth pulse, in which is encoded 
the protocol. A particular orientation of the register al- 
lows to suppress the dipole-dipole interaction between 
spins [HQj. 

The Hamiltonian of the system is 

1 n-1 j n-2 T 

h = --J2 - 2 E - E ^ (m) (*)e (ro) (*) 

i=o 2=0 m=l 

(1) 

with 

V (m \t) = — $>r exp{-i(^ + ^M )} + h . c0) 

1=0 

0( m )(i) equal to one during the mth pulse and zero oth- 
erwise, af' v ' z the usual Pauli operators for spin I and 
af~ = af ± iof . Due to the constant gradient of the per- 
manent magnetic field, the Larmor frequencies depend 
linearly on I, u>i — (I + l)a. By appropriately choosing 
the energy units we fix J = 1 throughout the paper so 
that the only relevant energy scales are and a. The 
basis states are chosen so that of \0)i = |0);. 

We will introduce the following notation for further 
discussion. Let the pulse P? c indicate a pulse with 
frequency vf c resonant with the flip of the i-th spin 
if its neighbors are in states "a" and "c". This will 
induce the resonant transition | . . . ai+xbiCi-x . . . ) — > 
| . . . cii+ibiCi-i . . . ), named T°" c , if the pulse is a 7r pulse 



(a,b,c € {0,1}). Note that for edge qubits, i.e. i € 
{0, n — 1}, only one superscript is needed. 

Operating IQC in the selective excitation regime, 
fi (m) < J < a, allows one to separate transitions in- 
duced by pulses into three sets: resonant, near-resonant 
and non-resonant according to the detuning A of the 
transition which is the difference between the frequency 
of the pulse and the energy difference of the states in- 
volved in the transition. If A is exactly equal to zero, 
the transition is called resonant (T°- c induced by the pulse 
P" c ), if A is of the order of J it is called near-resonant 
(Tf V induced by the pulse Pf c with {a',c'} ^ {a,c}), 
and if A is of the order of a it is called non-resonant 
(T" c induced by the pulse P°- c with i' ^ i). In the im- 
plementation of a protocol resonant transitions are the 
ones wanted, while near-resonant and non-resonant tran- 
sitions are a source of error. 

In the two level approximation a given unwanted 
transition with detuning A is induced with probability 

Q 2 n / 7T / aA 

P = fFTA^ sm ^2V 1+ f^J' (3) 

where p is a dimcnsionless duration of the pulse (for a 7r 
pulse p = 1 and for tt/2 pulse it is 1/2). The most prob- 
able transitions are the near-resonant ones and these can 
be suppressed as briefly described in the next paragraph. 

For Pj 10 (= P^ 1 ) pulses all near-resonant transitions 
have the same detuning so setting to A/V4A: 2 — 1 with 
k an integer suppresses these transitions. Since for near- 
resonant transitions A = O(J), Rabi frequency is for all 
pulses of the order of fl J/k. On the other hand, for 
Pf® and P/ 1 pulses near-resonant transitions have two 
different detunings therefore it is impossible to suppress 
both with a single pulse. This problem can be overcome 
adding an additional correcting P/ pulse. The combina- 
tion of these pulses in order to suppress all near-resonant 
transitions is called Q-pulse denoted by Qfj; when doing 
a p-K rotation of the ith qubit if neighbors are in states 
"a" and "c" . This method to eliminate near-resonant 
transitions is called generalized 2irk method. We refer 
the interested reader to Ref. |(| for further details, (im- 
pulses are the basic building blocks of gates, which in 
turn are the building blocks of algorithms such as QFT 
and IQFT. 

QFT for n = 4 qubits can be written as 

C^QFT = TA0B01B02B03A1B12B13A2B23A3. (4) 

There are in total n Hadamard A gates , n(n — l)/2 two- 
qubit B gates, B^-fe = diag{l, 1,1, exp (16 jk)}, with 9jk = 
tt/2 ^ and one transposition gate T which reverses the 
order of qubits (e.g. T|001) = 1 100)). In total there are 
n(n+l)/2 + l gates. IQFT algorithm Q for n = 4 qubits 
is given by 

£^iqft = TA0R01R02R03G01G02G03 

x A1R12R13G12G13A2R23G23A3, (5) 
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where Gy := RLBy. The R gate is defined by 
Rij\ . . . <u . . . bj . . .) := (-1)^ \...a i ...(cr i ®b j )...). In 
total there are n 2 + 1 gates in IQFT, i.e. roughly two 
times as many as for QFT. 

Recall that the implementation of quantum gates on 
IQC is easier in the interaction frame. Therefore, pulse 
sequences used in the paper implement the intended gates 
in the interaction frame. Each gate for QFT or IQFT 
(Eqs. (@J and (J5J) must in turn be implemented by sev- 
eral pulses (see appendix) . The number of pulses for QFT 
grows as ~ 18n 3 whereas it grows as ~ 54n 3 for IQFT. 
Note that this number can become very large, e.g. for 
IQFT and n = 10 one has 44541 pulses. Throughout 
the paper our basic unit of time will be either a gate (as 
written for instance in Eqs. an d ©) or a pulse. A 
single exception will be the paragraph discussing corre- 
lation function of intrinsic errors, where the basic unit is 
a Q-pulse, which is composed of one or two pulses. The 
reason is that Q-pulses are the smallest near-resonant 
corrected unit of generalized 2irk method. 

III. LINEAR RESPONSE THEORY 

As a criteria for stability we will use the fidelity F(t), 
defined as an overlap between a state ip(t) obtained by 
an evolution with an ideal algorithm and the perturbed 
ips(t) obtained by the perturbed evolution: 

F(t) = \{Mt)\m)\ 2 , (6) 

where \i/>(t)) = U(t)\^{0)) and \r/jg(t)) = U s (t)\i/)(0)). To 
simplify matters we will assume time t to be a discrete 
integer variable, denoting some basic time unit of an al- 
gorithm, like a gate or a pulse. The quantity measuring 
the success of the whole algorithm is the fidelity F(t) at 
t = T where T denotes the number of gates (pulses). 
One of the most useful approaches to studying fidelity is 
using linear response formalism in terms of correlation 
function of the perturbation, for a review see Ref. |l3j |. 
This approach has several advantages. First it rewrites 
the complicated quantity fidelity in terms of a simpler 
one, namely the correlation function, simplifying the un- 
derstanding of the fidelity. Second, the scaling of errors 
with the perturbation strength, Planck's constant and 
with the number of qubits is easily deduced. Further- 
more, as in practice one is usually interested in the regime 
of high fidelity, linear response is enough. 

First we will shortly repeat linear response formulas as 
they will be useful for our discussion later. Let us write 
an ideal algorithm up to gate t as U(t) 

U(t) = U t U t - l ...Ui, (7) 

where Ui is the i-th gate (pulse). If t = T we have a 
decomposition of a whole algorithm. 

The perturbed algorithm can be similarly decomposed 
into gates 

U s {t) = U 5 t U s t _ 1 ...Ul (8) 



Each perturbed gate Uj is now written as 

Uj = exp (-iSVj)Uj, (9) 

where Vj is the perturbation of j-th gate and S is a di- 
mensionless perturbation strength. For any perturbed 
gate one can find a perturbation generator V , such that 
the relation © will hold. Observe that the distinction 
into perturbation strength 5 and the perturbation gener- 
ator V in Eq. @ is somehow arbitrary. If one is given 
an ideal gate U and a perturbed one II s one is able to 
calculate only a product SV . This arbitrariness can al- 
ways be fixed by demanding for instance that the second 
moment of the perturbation V in a given state equals to 
1, {V 2 ) - {Vf = 1. 

To the lowest order in the perturbation strength fi- 
delity can be written as @ 

t 

F(t) = l-S 2 ]T C(h,t 2 ), (10) 

*1,*2 = 1 

where the correlation function of the perturbation is 

C(t u t 2 ) = (V tl (h)V t2 (t 2 )) - (V tl (h))(V t2 (h)) (11) 

with Vj(t) — U>(t)VjU(t) being the perturbation of j-th 
gate propagated by an ideal algorithm up to time t, i.e. 
in the Heisenberg picture. The brackets (•) denote the 
expectation value in the initial state. Throughout the 
paper we use random gaussian initial states and aver- 
age over many of them to reduce statistical fluctuations. 
Note that the time dependence of the correlation func- 
tion l|ll[l is due to two reasons: one is time dependence 
due to the Heisenberg picture (time index in brackets) 
and the second one is due to the time dependence of the 
perturbation itself (time index as a subscript), i.e. one 
has different perturbations Vj 7^ Vk for different gates 
j, k. The expression for the fidelity Eq. I|10|) is the main 
result of the linear response theory of the fidelity. From 
this one can see that decreasin g th e correlation sum (or 
even making it zero, see Ref. [lj|) W1U increase the fi- 
delity. In Ref. stability of QFT algorithm was consid- 
ered with respect to static GUE perturbation. Analyzing 
the correlation function they were able to design an im- 
proved QFT algorithm (IQFT) which increases fidelity. 

We are mainly interested in the fidelity F(T) at the 
end of an algorithm. The final time T in useful quantum 
algorithms depends on the number of qubits in a poly- 
nomial way, say as T oc n p . The power p depends on the 
algorithm considered and of course also on our decompo- 
sition of an algorithm into gates (pulses) . For QFT and 
IQFT algorithms with decomposition into gates, Eqs. © 
and (J5J, one has p = 2. On the other hand, for the imple- 
mentation of QFT on IQC one needs T oc n 3 (p — 3) basic 
electromagnetic pulses, as one is not able to directly per- 
form Hjk gates on distant qubits but has to instead use 
a number of pulses proportional to the distance between 
the qubits \j — k\. Now if the correlation function decays 
sufficiently fast, the fidelity will decay like F = 1 — S 2 an p 
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whereas in the case of slow correlation decay the fidelity 
will decay as F = 1 — S 2 cn 2p . In the extreme case of 
perturbations at different gates being statistically uncor- 
related (very fast decay of correlations) (VjVk) oc Sjk one 
obtains the exact formula F = exp(— S 2 n p ) In the 
limit of large quantum computer (large n) strongly cor- 
related static errors, giving slow decay of correlations, 
will therefore be dominant due to fast oc n 2p growth. 
When we will discuss errors caused by perturbations due 
to the coupling with the environment we will focus on 
static perturbations, meaning the same perturbation on 
all gates, Vk = Vj = V, as this component will dominate 
large n behavior. 



IV. ERRORS IN QFT 

Errors in an experimental implementation of QFT al- 
gorithm on an IQC can be of three kinds: (i) due to 
unwanted transitions caused by electromagnetic pulses 

(ii) due to coupling with external degrees of freedom and 

(iii) due to variation of system parameters in the course 
of algorithm execution. In the present paper we will dis- 
cuss only the first two errors. Errors due to electromag- 
netic pulses arc inherent to all algorithms on an IQC 
as we are presently unable to design pulse sequences for 
quantum gates without generating some unwanted tran- 
sitions albeit with small probabilities. This errors can 
be in principle decreased by going sufficiently deep into 
selective excitation regime but one must of course keep 
in mind the limitations of real experiments |l7j. Coupling 
with the "environmental" degrees of freedom is endemic 
in all implementations of quantum computers. As the 
environment will usually have many degrees of freedom 
we will model its influence on the quantum computer by 
some effective perturbation V c g given by a random ma- 
trix from a Gaussian unitary ensemble (GUE) Q. Note 
that coupling with the environment will generally cause 
non-unitary evolution of the central system. We expect 
quantum computation to be stable only on a time scale 
where evolution is approximately unitary, i.e. for times 
smaller than the non-unitarity time scale. Therefore we 
limit ourselves to unitary external perturbations. The 
third kind of errors due to the variation of system pa- 
rameters, like e.g. variation of Larmor frequencies due 
to the variation of the magnetic field is not considered in 
this paper. This does not mean they are not important. 
Let us consider a systematic error in the gradient of the 
magnetic field throughout the protocol (a — » a + Sa). De- 
manding that the error in the largest eigenphase at the 
end of the algorithm is much smaller than 1, one gets 
the condition a/Sa « an p+2 /ft. If one is in the selective 
excitation regime, this ratio can become very large and 
this puts a severe demand on the experiments. 

To ease up understanding we will first discuss intrinsic 
errors only, then we will discuss external errors only and 
finally we will combine both errors. 



A. Intrinsic Errors 

Let us consider probabilities of non and near-resonant 
transitions. For near-resonant and non-resonant transi- 
tions we have A 3> SI and the probability given by the 
perturbation theory is p oc (£!/A) 2 . In the 2wk method 
Rabi frequency Jl is given as 51 ~ J/k so the probabilities 
for near-resonant and non-resonant transitions are 

«■ « tub)) 2 ' < 12 » 

where p"° n denotes probability of a non-resonant tran- 
sition with A « a\j — 2 1 involving j-th and 2-th spin, 
one of which is a resonant one. The dependence of near 
and non-resonant errors on system parameters is there- 
fore different. 

For pulse sequences used to generate QFT or IQFT 
we always used a generalized 2irk method by which one 
can get rid of all near-resonant transitions. Therefore 
the only errors that remain are non-resonant ones. We 
first checked numerically that this is indeed the case by 
studying dependence of errors on system parameters by 
which one is able to distinguish near and non-resonant 
errors, Eq. 
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FIG. 1: Dependence of fidelity on physical parameters of the 
system. Empty points indicate variation of k with a = 100, 
filled points indicate variation of a with k = 128, both for 
n — 6. Full line is theoretical dependence of p non given by 
Eq. JTJ. 

As one can observe from Fig^the agreement with the 
theoretical p non (Eq. (|12fl ) is excellent thereby confirming 
that the only errors left are non-resonant ones. By using 
generalized 27rfc method we therefore decreased intrinsic 
errors by a factor of (a/ J) 2 as compared to ordinary 2irk 
method where there are still some near-resonant errors 
present. In order to have complete understanding of fi- 
delity decay due to intrinsic errors we have to understand 
scaling of these with the number of qubits. As we already 
discussed in section ITTT1 this depends on two things: how 
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strong the errors are correlated, giving possible scalings 
from n p to n 2p and on the increase of the perturbation 
strength with the number of qubits. Let us first dis- 
cuss the later. Under the assumption that the average 
transition probability (i.e. perturbation strength) for a 
non-resonant transition is the sum of all possible non- 
resonant transitions averaged over all possible resonant 
qubits, we can estimate 




with a some n independent constant. We can see that the 
perturbation strength does not grow with n asymptoti- 
cally, but the convergence to its limit is logarithmically 
slow. For small n the perturbation strength therefore 
will grow with n whereas it will saturate for large n. The 
second contribution to the n-dependence of the fidelity 
comes from the dynamical correlations between errors 
given by the correlation function l|llf> of the perturbation 
generator for non-resonant errors. We numerically cal- 
culated this correlation function in order to understand 
how the correlation sum and therefore fidelity, Eq. (|10|l . 
behaves as a function of n. 
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FIG. 2: (Color online) Correlation function for intrinsic errors 
in QFT for k = 128, a = 100, n = 4. The shading on time 
axes denotes the duration of different gates, Eq. J2J, and the 
time going from 1 to 543 runs over all Q-pulses. 

In Fig. [21 we show C(t\,t2) averaged over all Hilbert 
space. One can see that there are large 2-dimensional 
regions of high correlations in all parts of the picture. 
This means there are strong correlations between errors 
at different pulses and therefore the correlation sum will 
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FIG. 3: Fhe correlation sum of the same data as in Fig. [5] 
Note the decrease of the sum when the transposition gate is 
applied. The fidelity is in this linear response regime simply 
given by Eq. (1101 . Vertical lines indicate beginning of each 
gate (Eq. g}). 

likely grow as ~ n 6 as the number of pulses scales as n 3 
for our implementation of QFT. Similar results are ob- 
tained also for IQFT as can be inferred from Fig.0J One 
interesting thing to note is that during the application 
of the transposition gate at the end of the protocol the 
correlation sum starts to decrease at some point, nicely 
seen in Fig. and also visible in the correlation picture 
in Fig. as there are more negative than positive areas 
towards the end of the algorithm. This very interest- 
ing phenomena means that applying transposition at the 
end is advantageous (as compared to doing it classically 
for instance) as it will decrease non-resonant errors. We 
checked that this principle can not be exploited further 
by repeating transposition many times and by this de- 
creasing correlation sum even more. Still, this surprising 
behavior suggests that it might be possible to decrease 
non-resonant errors in a systematic way. 

To furthermore confirm predicted ~ n 6 growth of the 
correlation sum we calculated the dependence of intrinsic 
errors on n. This can be seen in Fig. 0] where we plot 
1 — F(T) as a function of n for QFT and IQFT and for 
two different sets of parameters, one for k = 128, a = 100 
giving large errors and one for k = 1024, a = 1000. One 
can see that asymptotically for large n the dependence is 
indeed n 6 but the convergence to this behavior is fairly 
slow, one needs of the order of n = 7 or more qubits. 
This slow convergence we believe is due to the logarithmic 
convergence of the perturbation strength (Eq. (|13fl b To 
get exact coefficients in front of n 6 dependence we fitted 
dependences of errors in Fig. 01 with a polynomial in n 
using at most two nonzero terms. Defining polynomials 
in the linear response regime as s ln = (1 — F)(ka/ J) 2 one 
gets for QFT and IQFT 

SqftH = 280n 6 - 660n 5 
s jn FT (n) = 1300n 6 - 2100n 5 . (14) 

Both expressions are good for n > 5 and superscript "in" 
denotes intrinsic errors. Beyond the linear response the 
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FIG. 4: Dependence of fidelity on the number of qubits. 
Empty symbols indicate data for k = 128 and a — 100 while 
filled symbols are for k — 1028 and a = 1000. Circles indi- 
cate QFT and squares IQFT. In the presence of only intrinsic 
errors, IQFT does not improve fidelity. Full lines show asymp- 
totic n 6 dependence. 



exponential dependence is frequently justified and one 
has 



F = exp 



J 

ka 



s in {n) 



(15) 



Large coefficients of polynomials in Eqs. (|14fl are due 
to large number of pulses. The maximum possible de- 
pendence in the case of no decay of correlation function 
(see discussion at the end of section could be T 2 
and therefore the leading terms in polynomials l|14|l ex- 
pressed by the number of pulses are Sq FT ~ 0.8Tq FT and 

s 1qft ~ 0.5Tiq FT . Therefore relative to the number of 
pulses IQFT slightly decreases non-resonant errors but in 
the absolute sense QFT is better simply because it has 
only one third as many pulses as IQFT and the coefficient 
in front of n 6 (Eq. (114(1 ) is thereby smaller. If only intrin- 
sic errors in the generalized 2irk method are concerned 
QFT is always more stable than IQFT. Note that the 
intrinsic errors due to non-resonant transitions for QFT 
grow as ~ T 2 (~ n 6 ) whereas in previously studied "sim- 
ple" algorithms, for instance entanglement protocol Q, 
they grow only as the first power of the number of gates 
~ T. This means that QFT is much more sensitive to 
intrinsic errors. 



B. External Errors 

In order to study only external errors we set through- 
out this section parameters to k = 1024 and a = 1000, 
for which intrinsic errors are much smaller than external 
ones. 

External error will be modeled by the perturbation V 
(Eq. @) chosen to be a random hermitian matrix from 
a GUE ensemble. To facilitate comparison with previous 



results on IQFT || we will make perturbation after each 
quantum gate, except for the last transposition gate T, 
Eq. after which we do not make perturbation. So for 
QFT we make n(n + l)/2 perturbations, while for IQFT 
we make n 2 perturbations. One other possible choice 
would be to make perturbations after each pulse. We 
will discuss this possibility at the end of this section. For 
now let us just say that qualitatively the results are the 
same as if doing perturbation after each gate, one just 
has to rescale perturbation strength like <5 ga t c oc nS pu \ sc 
as there arc effectively 0(n) perturbations (pulses) per 
gate. 

The implementation of QFT on IQC is written in the 
interaction picture. As the static perturbation is the 
worst, meaning it will asymptotically in large n limit be 
dominant, we will concentrate only on static perturba- 
tion, i.e. the same perturbation for all gates (pulses) 
Vj = Vk = V. There are still two possibilities, cither 
making static perturbation in the interaction frame or 
making it static in the laboratory frame. Let us first 
discuss the later case. If we make static perturbation 
in the laboratory frame, we can of course transform it 
to the interaction frame by a unitary transformation 
W(t) = exp(— iHot) given by the time independent part 
Hq of the Hamiltonian Eq. QJ. This transformation 

exp (-WKtot(t)) := W\t) exp (-iSV la b)W(t), (16) 

will result in the perturbation in the interaction frame 
Vi n t(i) being time dependent. As the transformation 
to the laboratory frame W(t) in the selective excitation 
regime involves large phases, perturbations at different 
gates will tend to be uncorrelated due to averaging out of 
widely oscillating factors in the correlation function i|llll . 
Therefore in the first approximation one can assume 
C{t\ 1 t2) = <5t!,t 2 and the fidelity will in this extreme case 
of uncorrelated errors decay as F(T) = exp (— 8 2 T) 
On the other hand, if we make perturbation to be static 
in the interaction frame and the correlation function is 
not Kronecker delta in time, the fidelity will decay faster 
with the number of qubits (or T). To numerically con- 
firm these arguments, we show in Fig. [S] dependence of 
fidelity with the number of qubits n for the two cases dis- 
cussed, static perturbation in the interaction frame and 
static perturbation in the laboratory frame. Polynomial 
fitting of n dependence for QFT gives 



sf a u b c (n) = 0.47n 2 + 1.41n- 2.42 
s f„ U t°H = 0.45n 3 - 0.42n 2 + 0.58?i. 



(17) 



The fidelity due to external GUE errors is given as 

F = exp(-<5 2 s suo (n)), (18) 



with the appropriate s gue (n) from Eq. (|17fl . Note that 
s?" t c (n) grows faster than sf^(n) as argued. Observe 
also that for the static perturbation in the laboratory 
frame using the assumption of uncorrelated errors in the 
interaction frame we predicted sf^ w T w n 2 /2 for QFT 
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which is remarkably close to the numerically observed 
value Eq. {T7JI- For IQFT and the application of GUE 
perturbation in the laboratory frame one gets a similar 
result with the leading term sf^(n) ~ 1.12rt 2 . One can 
write an arbitrary time dependent perturbation in the 
interaction frame as a Fourier series and for large n the 
static component will always prevail. Therefore, from 
now on we will exclusively discuss only static perturba- 
tions in the interaction frame. 




3 4 5 6789 10 



n 

FIG. 5: Dependence of 1 — F in QFT algorithm on the num- 
ber of qubits for static perturbation in the interaction frame 
(empty points) and in the laboratory frame (full points), with 
S = 0.04. Lines are best fitting polynomials (see Eq. HTH h 

The dependence of errors due to GUE perturbation in 
the case of QFT and IQFT has already been derived [lg| . 
For IQFT numerical fitting in our case gives dependence 

sfq FT (n) = l.Sln 2 + OMn - 3.73. (19) 

Here the perturbation is again static in the interaction 
frame as is the case throughout the paper with a single 
exception being the previous paragraph. Dependence of 
fidelity in both cases for QFT and IQFT can be seen in 
Fig- EH together with the theory (Eq. iJTHJ) using polyno- 
mials iJTJJJ and ifPTjl, 

Observe that IQFT for n > n cr it = 3 is better than 
QFT despite having more gates and therefore applying 
perturbation on it more times (for n < 3 QFT is slightly 
better). What is important is that the dependence of 
errors on n is also different, ~ n 3 for QFT, but only ~ n 2 
for IQFT. This means that asymptotically IQFT is much 
more stable against GUE perturbations than ordinary 
QFT. 

Finally, let us discuss what happens if we make static 
GUE perturbations in the interaction frame after each 
pulse, and not after each gate as done so far. The 
product of two operators exp (-iSV)U can be written as 
exp(-i5V)U = Uexp(-iSV(l)), with V(l) := WVU. 
Using this expression for all errors within a single gate 
and bringing them all to the beginning of the gate we get 

exp {-iSV)U r • • ■ exp {-\8V)Ui w 
Ur-'-Ut exp {-iS[V(l) + ■■■ + V(r)\), (20) 
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n 

FIG. 6: Dependence of fidelity on the number of qubits for 
QFT (empty symbols) and IQFT (filled symbols) algorithms 
(8 = 0.04). Curves are theoretical prediction Eq. I|18jl using 
polynomials from Eqs. 1191 and I17H . 

where r is number of pulses constituting a gate. This 
means that the application of the perturbation after each 
pulse is to the lowest order in S equivalent to the appli- 
cation of the effective perturbation SY^jV(J) after the 
gate. Of course now the perturbation is explicitly time 
dependent. But for a GUE matrix acting on a whole 
Hilbert space individual pulses will do transformations 
on an exponentially small subspace (i.e. on one qubit) 
and therefore one might expect that effectively one can 
write S^2j V(J) ~ rV e ff, where V e s is some effective ran- 
dom matrix independent of the gate. As in our case 
of QFT on IQC we have on average oc n pulses for a 
single gate we can predict that doing perturbation with 
strength <5 pu i so after each pulse is approximately equal 
as doing perturbation of strength (5 ga t ~ ^^puise after 
each gate. In order to confirm these expectations we did 
numerical experiments with the results shown in Fig. 
Fitting polynomial in the dependence of fidelity, Eq. iJTSJ, 
for QFT and IQFT gives in this case 

s| U p T (n) = 4.86n 5 + 35.8n 4 
s^p T (n) = 25.6n 4 + 606n 3 . (21) 

The leading dependence of n 5 for QFT and n 4 for IQFT 
nicely agrees with our rescaling prediction S ga t c ~ n<5 pu i SG . 
IQFT is asymptotically again better than QFT as the er- 
rors grow slower with the number of qubits. The crossing 
point between the two in this case happens at n crit = 10, 
whereas in the case of perturbation after each gate we 
had n cr it = 3. This confirms that doing GUE perturba- 
tion after each pulse is qualitatively the same as doing 
it after each gate, only the crossing point between QFT 
and IQFT changes and of course also the dependence of 
errors on n changes, simply due to the different number 
of applied perturbations. If perturbation strength 5 is 
properly rescaled, the n dependence is the same in both 
cases. 

Up to now we discussed intrinsic errors and external 
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FIG. 7: Dependence of F on the number of qubits for the 
static GUE perturbation after each pulse with S = 5 ■ 10~ 4 . 
Empty symbols are for QFT and filled symbols are for IQFT. 
Curves are theoretical prediction Eq. 1181 using best fitting 
polynomials given by Eq. 112 II . 



FIG. 8: Fidelity for QFT (empty symbols) and IQFT (filled 
symbols) algorithm and GUE perturbation after each gate. 
System parameters are k — a = 200 and S = 0.04 (intrinsic 
and external errors are comparable in size). Full curves are 
theoretical predictions for F given by Eq. 1221 . 



errors separately. The next question is of course, what 
happens if both errors are present at the same time and 
are of similar strength? 



C. Intrinsic and External Errors combined 

If both kinds of errors are present, a first naive guess 
would be that they just add, 





' J " 


- 


ka 



cxp - f- s in (n)-5 2 s SUB (n) , 



^iboth _ pinjpgue _ 



with the appropriate polynomials s ln (n) and s gue (n) 
given in previous Eqs. I|14|). (|17fl and (|19fl . In the linear 
response regime this formula means that both errors are 
uncorrelated, i.e. their cross-correlations are zero. This 
is easy to proof using properties of GUE matrices. Let us 
calculate cross-correlation function between V m (ii) and 
V sue (t2) averaged over GUE ensemble. Written explic- 
itly one has to average products of the form V™Vp? c ', 
where V guo is a GUE matrix. As this expression is linear 
in V guc it averages to zero, (V^V^ ) guc = 0, thereby 
explicitly confirming a simple additions of both errors. 
Of course in real experiments we are not averaging over 
GUE ensemble but are taking one definite representative 
member of it. But for large Hilbert space the expecta- 
tion value of a typical random state and one particular 
GUE matrix is "self-averaging" and will be equal to the 
ensemble average. 

Let us check the theoretical prediction for fidelity 
Eq. Q22JI with a numerical experiment. We again apply 
GUE perturbation after each gate. The results together 
with the theoretical prediction Eq. Il22[l are in Fig. [S] 
The agreement between the theory and the experiment 
is good also beyond the linear response regime. Please 
note that we deliberately choose parameters so that both 



QFT and IQFT give similar fidelity in order to also see 
the crossing of the two curves within the shown range of 
n. Given fixed S and ka, QFT is always better for large 
n because intrinsic errors will prevail over external ones, 
due to their fast n 6 growth. But still, for intermediate 
n's IQFT can be better that QFT as seen in Fig. 

Now we are equipped with understanding of errors in 
QFT and IQFT due to external GUE perturbation and 
intrinsic errors so we can make some predictions regard- 
ing ranges of experimental parameters ka, 5, n for which 
the fidelity will be high enough. Interesting question for 
instance is, when is IQFT better than QFT? To find that, 
we set -Fqft = -Fiqft with F's given by Eq. (|2*2l . This 
results in the condition 



^crit — -, — -I 

ka 





„in 


S QFT 


IQFT 


„guo 





'QFT IQFT 



(23) 



For S > S CTlt IQFT is better than QFT. In Fig. we 
show curves of constant fidelity for n — 5. They are com- 
posed of two parts, above the line for 5 CI [t IQFT is better 
than QFT, and below vice versa. Two characteristic fea- 
tures are also vertical and horizontal asymptotes of the 
curves of constant fidelity. The vertical asymptote means 
that for fixed n, even if 6 = 0, we must have ka larger 
than some critical value determined just by intrinsic er- 
rors, in order to have given fidelity. Horizontal asymptote 
for high ka means that if 5 is larger than some critical 
value, increasing ka will not help to improve fidelity. In 
Figs. EI| an d ^2 we show similar plots, only these time 
one of the axes is dependence on n. For instance, from 
Fig. ^3 on can see that having ka = 10 5 , the maximum 
number of qubits is n ss 12 if we want to have fidelity 
larger than 0.9 (even if S = 0). This unfavorable growth 
of required ka oc n 3 in order to have a fixed fidelity is 
due to ~ n 6 growth of intrinsic errors. It would therefore 
be advantageous to find a way to suppress errors due to 
non-resonant transitions |l6j . 
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FIG. 9: (Color online) Dependence of fidelity on system pa- 
rameter ka and GUE perturbation strength 8 for n = 5. Full 
curves of constant fidelity are composed of two parts corre- 
sponding to QFT or IQFT. Above the thick line for 5 C rit IQFT 
is better and below QFT is better. Dotted curves of constant 
fidelity below this line are for IQFT and dashed lines above 
are for QFT. The shaded region corresponds to the region of 
fidelity larger than 0.9. The plus symbol shows the position 
of parameters for Fig. |H| 
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FIG. 11: (Color online) Fidelity dependence on ka and n for 
a fixed 8 = 0.01. For the explanation of various curves see 
the caption to Fig. [§| 



algorithms (having 0(n) gates), where the growth is lin- 
ear in time. The immediate question is whether this be- 
havior is general for algorithms having more than 0{n) 
gates. This very fast growth with n is a consequence 
of strong correlations between errors at different pulses 
and puts a severe demand on experimental requirements. 
Therefore it would certainly be desirable to find a way 
to suppress also non-resonant errors. We also consider 
perturbations due to coupling with external degrees of 
freedom modeled by a random GUE matrix. To sup- 
press this kind of errors we show that it is advantageous 
to use an improved QFT algorithm, for which the errors 
grow only as ~ n 2 , whereas they grow as ~ n 3 for or- 
dinary QFT. By a combination of both techniques, the 
generalized 2irk method and improved QFT algorithm, 
we are able to make implementation of QFT stable in a 
much wider range of parameters. 



FIG. 10: (Color online) Fidelity dependence on 8 and number 
of qubits n for a fixed value of ka = 10 5 . For the explanation 
of various curves see the caption to Fig. |3] 



V. CONCLUSIONS 

We analyzed two possible errors in the implementa- 
tion of QFT on IQC working in the selective excitation 
regime. We consider: (i) intrinsic errors due to unwanted 
transitions caused by pulses, (ii) external errors due to 
coupling with the external degrees of freedom. We care- 
fully analyze their dependence on system parameters and 
on the number of qubits. To diminish intrinsic errors we 
use the generalized 2nk method by which we are able to 
suppress all near-resonant transitions, with only much 
smaller non-resonant transitions remaining. We then 
study these non-resonant errors in QFT algorithm and by 
using correlation function formalism explain their growth 
with time as ~ T 2 , in contrast to so far studied "simple" 
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APPENDIX: QFT AND IQFT 
IMPLEMENTATION ON THE ISING QUANTUM 
COMPUTER 

To implement the protocol with high fidelity we use 
Qfp pulses derived in Ref . @ , which completely suppress 
all near-resonant errors. Phases of Q-pulses composing 
a gate must be chosen correctly so that the gate works 
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on an arbitrary state. The protocols implementing CNy 
(control not gate) and Nj (not gate) can be found in 
sections 7.1-7.3 of Ref. Q. 

In order to complete QFT and IQFT we still need to 
implement the R T , R, A, B and T gates. We can decom- 
pose R, R* and T gates into simpler pieces: 



R 



NjCNj,-NjZ,-, 



K,', = NiZjCN^-Ni, and 

q q-i 

T = 1 I 1 I S / V 1 

i=lj=l 



(A.l) 
(A.2) 

(A.3) 



with Sjj = CNijCNjiCNy the swap gate, Z = 
diag{l,— 1} the a z gate and each term in the product 
in Eq. (|A. 3|) is placed at the left of the sub-product (e.g. 
lli=o Di — D2D1D0). Therefore, the only gates left to 
design are A, B and Z. 

The phases of Q-pulses can be expressed in terms of 
angles 9 P , a p , 9 P , (3 P and j p |6| which are given by 



e p = -K^ki - pV4, (a.4) 

«p = | V /^+V/4 ! (A.5) 

(A.6) 
(A.7) 
(A.8) 



tan0 o = — - 



2a r , 



tana. 



tan/3 p = — - — tana p cos9 p , 

LOi p 

lp = ^( 7r fc p )2_( 7r + /3p )2 



We use notation of angles without subscripts denoting 
angles for tt pulses i.e. 9 = 9\ and set fci/2 = 2k. 
The Hadamard gate can now be expressed as 



A 



, - Qf{9i)Qf{v>2)Qf{^)QfM 



(A.9) 



QT^)Q^/2), 
for intermediate qubits and 

Aj = Q°(^ 6 )Q> 7 )Q°i {v&)Q% (tt/2), (A.10) 
for edge qubits, with 



pi = -2 (f + 7| +6 

¥>3 = -2 (0 + 7 - 71 - 
055 = tt/2 — 6>i — 29 1 , 

2 2 

V? 7 = -0 + 0±, 



¥>2 


= -0- 


26, 




= tt/2- 


-2 7 


<P6 


= -0- 


d b 


¥>8 


= n/2- 


- 29 1 



For neighboring qubits (\i — j'| = 1) the B gate can be 
written as, 



Bii = Q I 11 (0)Q I 10 (0)g° (0)g} (0)Q)°(^ 1 )Q° (0) 

Q) (v>dQ?(W?(vs), (A.12) 

for intermediate qubits and for edge qubits (i or j G 
{0,n- 1}) it is 



B, 



Qj(0)Q?(0)Qf(0)Qf(0)Q?°(0) 

Qf(w)Qite7)Q?(w0Qf(o) 



(A.13) 



QffoOQf (0)Qffo> 10 ). 
Angles for B gates are 



l Pi = -27 - 361 + 29, y> 2 = cj)/2 - 2 7 - 66 1 , 

<^3 = 0/4-7r/2, ^4 = -(fx, 

<p 5 = 0/2 + 2 7 + 66>, ^ 6 = 0/2-67- 120 + 49, 

959 = -2i/3i, (^10 = 0/2-27 + 49, 

(A.14) 

and — 7r /2- For distant qubits (\i — j\ > 1) it is nec- 
essary to use swap gates to bring z-th and j-th qubits to 
neighboring positions, then apply B protocol for neigh- 
bor qubits and finally take them back to their original 
positions using swap gates. The angle in Eq. (|A.14|I is 
in this case = 7r/2' J " — Finally the Z gate is expressed 
as 



Zj = Q} 1 (0)Q) (0)QO°(0)Q} 1 (7r/2)Q} (7r/2)Q™( 7 r/2) 



for intermediate qubits and 

Zj = Q}(0)Q°(0)Q}(V2)Q J °(7r/2) 



(A.15) 



(A.16) 



for edge qubits. Counting the number of all pulses for 
QFT and IQFT one gets 



Tqft = 18n 3 - 16n 2 - 49n + 57 (A.17) 
T IQFT = 54n 3 - 86n 2 - 105n + 191. (A. 18) 



(All) 
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